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Abstract 

Measured responses from visual cortical neurons show that spike times tend to be 
correlated rather than exactly Poisson distributed. Fano factors vary and are usually 
greater than 1 due to the tendency of spikes being clustered into bursts. We show 
that this behavior emerges naturally in a balanced cortical network model with ran- 
dom connectivity and conductance-based synapses. We employ mean field theory 
with correctly colored noise to describe temporal correlations in the neuronal activ- 
ity. Our results illuminate the connection between two independent experimental 
findings: high conductance states of cortical neurons in their natural environment, 
and variable non-Poissonian spike statistics with Fano factors greater than 1. 
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1 Introduction 



Neurons in primary visual cortex show a large increase in input conductance 
during visual activation: in vivo recordings (see, e.g., [1]) show that the con- 
ductance can rise to more than three times that of the resting state. Such high 
conductance states lead to faster neuronal dynamics than would be expected 
from the value of the passive membrane time constant, as pointed out by Shel- 
ley et al. [2]. We use mean field theory to study the firing statistics of a model 
network with balanced excitation and inhibition and observe consistently such 
high conductance states during stimulation. 

In our study, we classify the irregularity of firing with the Fano factor F, 
defined as the ratio of the variance of the spike count to its mean. For tem- 
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porally Tincorrclated spike trains (i.e., Poisson processes) F = 1, while F > 1 
indicates a tendency for spike clustering (bursts), and F < 1 points to more 
regular firing with well separated spikes. Observed Fano factors for spike trains 
of primary cortical neurons during stimulation are usually greater than 1 and 
vary within an entire order of magnitude (see, e.g., [3]). We find the same 
dynamics in our model and are able to pin-point the relevant mechanisms: 
synaptic filtering leads to spike clustering in states of high conductance (thus 
F > 1), and Fano factors depend sensitively on variations in both threshold 
and synaptic time constants. 



2 The Model 



We investigate a cortical network model that exhibits self-consistently bal- 
anced excitation and inhibition. The model consists of two populations of 
neurons, an excitatory and an inhibitory one, with dilute random connectiv- 
ity. The model neurons are governed by leaky integrate-and-fire subthreshold 
dynamics with conductance-based synapses. The membrane potential of neu- 
ron i in population a (a = 1,2 for excitatory and inhibitory, respectively) 
obeys 

^ = -9L<{t) -EIlg'MKit) - H). (1) 

b=0 j=l 



The first sum runs over all populations b, including the excitatory input pop- 
ulation representing input from the LGN and indexed by 0. The second sum 
runs over all neurons j in population b of size A*";,. The reversal potential V5 
for the excitatory inputs (6 = 0, 1) is higher than the firing threshold, the one 
for the inhibitory inputs (V2) is below the reset value. The constant leakage 
conductance ql is the inverse of the membrane time constant r^. 

The time dependent conductance Qabit) from neuron j in population b to 
neuron i in population a is taken as 

gl^t) = J^J2^^p{-{t-ti)/n)Q{t-ti) (2) 



if there is a connection between those two neurons, otherwise zero. The sum 

runs over all spikes s emitted by neuron j, is the synaptic time constant for 
the synapse of type b (excitatory or inhibitory), and is the Heavyside step 
function. denotes the average number of presynaptic neurons in population 
b. We followed van Vreeswijk and Sompolinsky [4] by scaling the conductances 
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with Xj \/K\) so that their fluctuations are of order one, independent of network 
size. 



3 Mean Field Theory 

We use mean field theory to reduce the full network problem to two neurons: 
one for each population. This method is exact in the limit of large popu- 
lations with homogeneous connection probabilities [5]. The neurons receive 
self-consistent inputs from their cortical environment, exploiting the fact that 
all neurons within a population exhibit the same firing statistics due to ho- 
mogeneity. The time dependent conductance described in (2) can then be 
replaced by a realization of a Gaussian distributed random variable gah with 
mean 

(^ab) = glb^K^bn, (3) 

and covariance 

{5gab{t) 5ga,{t')) = {gl,f{l - K,/N,) C,{t - t'), (4) 

Here, is the firing rate of the presynaptic neuron b, and Cb{t — t') is the 
autocorrelation function of its spike train. A simple approximation of the 
autocorrelation, hke the one used by [6] and [7], is to assume gab{t) to be 
temporally uncorrelated (i.e., white noise), in which case it simplifies to Cb{t — 
t') — Tb S{t — t'). The term (1 — Kb/Nb) is a correction for the finite connection 
concentration Kb/Nb and can be derived using the methods of [8]. 

The self-consistent balance condition is obtained by setting the net current 
in (1) to zero when the membrane potential is at threshold 6a and the con- 
ductances have their mean values (3). In the large Xb-limit, it reads 

jZglb\fKbrb{ea-Vb)^Q. (5) 

6=0 

The distribution of the variables gab can be calculated numerically using an 
iterative approach [9]. One starts with a guess based on the balance equa- 
tion (5) for the means and covariances and generates a large sample of specific 
realizations of gab{t), which are used to integrate (1) to generate a large sample 
of spike trains. The latter can then be used to calculate new estimates of the 
means and covariances by applying (3) and (4) and correction of the initial 
guess towards the new values. These steps are repeated until convergence. 
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Fig. 1. Autocorrelation functions for white noise (blue) and colored noise (red). The 
white noise approximation underestimates the amount of temporal correlation in 
the neuron's firing. 

4 Results 

For the above described model, we chose parameters corresponding to popula- 
tion sizes of 16,000 excitatory neurons and 4,000 inhibitory neurons, represent- 
ing a small patch of layer IV cat visual cortex. The neurons were connected 
randomly, with 10% connection probability between any two neurons. The 
firing threshold was fixed to 1, excitatory and inhibitory reversal potentials 
were set to +14/3 and —2/3, respectively, and the membrane time constant 
Tm = gj} was 10 ms. For the results presented here, the integration time step 
was 0.5 ms. 

Figure 1 illustrates the importance of coloring the noise produced by intra- 
cortical activity. The white noise approximation underestimates both the cor- 
relation times and the strength of the correlations in the neuron's firing: its 
autocorrelation (blue) is both narrower and weaker than the one for colored 
noise (red). 

Fano factors vary systematically with both the distance between reset and 
threshold and the synaptic time constant r^. Non-zero synaptic time constants 
produced consistently Fano factors greater than one. We varied the reset be- 
tween 0.8 and 0.94 and between and 6 ms, which resulted in values for F 
that span an entire order of magnitude, from slightly above 1 to approximately 
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Fig. 2. Fano factors for a range of reset values and synaptic time constants Tg- 
Longer synaptic time constants lead to increased clustering (bursts) of spikes, which 
is reflected in higher Fano factors. 

10 for Tg > 2 ms (see Figure 2). 



5 Discussion 



In all our simulations, we observed that the membrane potential changed on a 
considerably faster time scale than the membrane time constant = 10 ms. 
This behavior is only observed if conductance-based synapses are included in 
the integrate-and-fire neuron model. To understand this phenomenon, it is 
convenient to follow the notation of Shelley et al. [2] to rewrite the equation 
for the membrane potential dynamics (1) in the following form: 

^ = -g^(t)Mt)-Vsit)), (6) 



with the total conductance grit) = Ql + J2b9ab(t), and the effective reversal 
potential Vs{t) = gT{t)^^J2b9ab{t)Vb. The membrane potential Ua{t) follows 
the effective reversal potential with the input dependent effective membrane 
time constant gT{t)~^ . The effective reversal potential changes on the time 
scale of the synaptic time constants, which are up to five times shorter than 
in our simulations. However, if the effective membrane time constant is shorter 
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Fig. 3. The membrane potential u{t) (red) follows the effective reversal potential 
Vs{t) (blue) closely, except for detours when the neuron is reset due to firing. The 
membrane potential recovers fast enough to spike several times while Vs{t) stays 
above threshold, thus producing bursts of spikes. Here, the threshold is set to 1 and 
the reset to 0.94. 

than the synaptic time constant due to a large enough total conductance, then 
Uait) can follow Vs{t) closely, as observed in our simulations (see Figure 3). 

In high conductance states, the firing statistics are strongly influenced by 
synaptic dynamics (see Figure 2). This is in contrast with strictly current 
based models, where the neuron reacts too slow to reflect fast synaptic dy- 
namics in its firing. The 'synaptic filtering' of arriving spikes leads to temporal 
correlations in Vs{t) and thus to temporal correlations (by way of spike clus- 
tering) in firing. Therefore, the model neurons receive temporally correlated 
input rather than white noise. For this reason, in mean field models dealing 
with conductance based dynamics, coloring the noise is important to arrive 
at the full amount of temporal correlation in firing statistics (see Figure 1). 
We confirmed these considerations by running simulations without synaptic 
filtering (r^ = 0). As expected, intra-cortical activity became uncorrelated and 
the white noise approximation produced the same result as coloring the noise 
correctly. In that case, Fano factors stayed close to 1 (see Figure 2), i.e, no 
tendency of spike clustering was observed. 

Previous investigations showed that varying the distance between threshold 
and reset in balanced integrate-and-fire networks has a strong effect on the 
irregularity of the firing [10] . By including a conductance-based description of 
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synapses, wc were now able to show the importance of synaptic time constants 
on firing statistics, even if they are several times smaller than the passive 
membrane time constant: Synaptic filtering facilitates clustering of spikes in 
states of high conductance. 
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